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Abstract 

In this paper, we study a free boundary problem for compressible spherically symmetric 
Navier-Stokes equations without a solid core. Under certain assumptions imposed on the 
initial data, we obtain the global existence and uniqueness of the weak solution, give some 
uniform bounds (with respect to time) of the solution and show that it converges to a sta- 
tionary one as time tends to infinity. Moreover, we obtain the stabilization rate estimates of 
exponential type in L°°-norm and weighted i? 1 -norm of the solution by constructing some 
Lyapunov functionals. The results show that such system is stable under small perturbations, 
and could be applied to the astrophysics. 

Keywords: Compressible Navier-Stokes equations; density-dependent viscosity; free bound- 
ary; existence; uniqueness; asymptotic behavior 

1 Introduction. 

We consider the compressible Navier-Stokes equations with density-dependent viscosity in R n (n > 
2), which can be written in Eulerian coordinates as 



d T p + V • (pu) = 0, 

d T (pu) + V • (pu gu)| VP = div(p(Vu + Vu ] )) + V(Adivu) - pf \ 



■T^ , „ MJI ^ .7 (1-1) 



Here p, P, u = (ui,...,u n ) and / are the density, pressure, velocity and the external force, 
respectively; p = p(p) and A = \(p) are two viscosity coefficients. 
In this paper, the initial conditions are 

p(C0) = p o (r), re [0,6], (1.2) 

u(£,0) = u (r)^,r G (0,6], u(£,0)| f=0 = «o(0) = 0, (1.3) 
where r = |£| = + • • • + and 6 > is a constant, the boundary condition is 

Up- Adivu)Id - /i(VS + Vn T )} • n = P r n, £ e dQ T (1.4) 
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where d£l T = ip(dQo,T) is a free boundary, n is the unit outward normal vector of dVt T and 
Pp > is a external pressure. Here, 8Qq = {£ G M n : |£| = 6} is the initial boundary and ?/> is 
the flow of u: 

f ^(f.r) = n(V(£r),r), £ G M n , f 5) 

\v>(£o) = e 

To simplify the presentation, we only consider the famous polytropic model, i.e. P(p) = Ap 1 
with 7 > 1 and A > being constants. And we assume that the viscosity coefficients \i and A 
are proportional to p e , i.e. p(p) = c\p and X(p) = C2p e where c\,C2 and are three constants. 

For the initial-boundary value problem (jl.l|) - (jl.4p . we are looking for a spherically symmetric 
solution (p, u) : 

PiC, r) = p(r, r), u(£, r) = u(r, r)-, 

r 

with the spherically symmetric external force 

f = f(m,r,r)-, m(p,r) = / p(s,T)s n ~ 1 ds, r>0 

r Jo 

and dfl T = {£g M n : |£| = 6(t),&(0) = 6,^(r) = w(6(r),r)}. 
Then (p, u) (r, r) is determined by 

f 9 T p + 3 r (pu) + ^V = 0, 

p(<9 T u + ud r u) + d r P = (A + 2 / u)(a2 r u + ^-d r u - ^±u) (1.6) 

+2d r pd r u + d r X(d r u + — p/, 

where (r, r) E (0, 6(r)) x (0, oo), with the initial data 

(p,n)| T=0 = (p ,u )(r), <r < b, (1.7) 

the fixed boundary condition 

U\ r =0 = 0, (1.8) 

and the free boundary condition 

P -2pd r u- \(d r u + — -u] 1 = P r , (1.9) 

V r / J r=6(r) 

where 6(0) = 6, b'{r) = n(6(r),r). 

Additionally, we assume the external force f(m,r, r) and external pressure Pt(t) G C 1 (M+) 
satisfy 

i^(T) = P 0O + AP(r), /(m,r,T)=/ 00 (m,r) + A/(m,r ) T), (1.10) 
for all r > and r > 0, with 

/oo(m,r) = —j, m(p,r) = / p/ 1 " 1 ^, A/(m,r,r) G C 1 ^ x R+ x E + ) (1.11) 
r Jo 

l|A/(-,-,r)|| LO o (R+xR+) <A(r), ||(a r A/,a T A/)(-,.,r)|| Loo(M+xR+) </ 2 (r), (1.12) 

/!, AP € L™ n L 1 ^), (AP)', / 2 G L 2 (M+), (1.13) 

where M+ = [0,oo), Px, and G are two positive constants, perturbations (AP, A/) tend to as 
t — > oo in some weak sense, /oo is the precise expression for its own gravitational force and Af 
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expresses the influence of the outside gravitational force, in the astrophysical case (with spherical 
symmetry). Pp also could express the influence of the surface tension force on the free boundary. 
This system can be treated as a simple model of one fluid in f2 T , whose evolution is influenced 
by the gravitational force and the external pressure generated by the other substance in M n \f2 T . 
We study the stabilization problem of such system, which could be applied to the astrophysics. 
Now, we consider the stationary problem, namely 

(P(Poc))r = -Poo/oo(™.(/Ooo,'r),r-) (1.14) 
in an interval r £ (0, 1^) with the end satisfying 

P(Poo(loo)) = Poo, (1-15) 
loo rb 

Poo r n ~ l dr = M := / p Q r n ~ l dr. (1.16) 



>o Jo 
The unknown quantities are the stationary density p^ > and free boundary > 0. If 

In — 2 2-n 2 

7 = and Gn~ M» < 2A (1.17) 

n 

or 

2n - 2 , 

7> , (1.18) 

n 

from Proposition 12.51 we know that there exists a unique solution (poo^oo) to the stationary 
system fjl . 14|) - ()1 . 16[) . satisfying < p < Poo(0 < p < oo, (poo)r(O < 0, < r < with 

loo < +00. 

To handle the free boundary problem (jl.6p - (jl.9p . it is convenient to reduce the problem in 
Eulerian coordinates (r, r) to the problem in Lagrangian coordinates (x,t), via the transforma- 
tion: 

x= f ' y n - l P {y,T)dy, t = r. (1.19) 



o 



Then the fixed boundary r = and the free boundary r = b{r) become 

Mr) rb 

x = andx= / y n l p(y,T)dy= / y n l p (y)dy = M, 
Jo Jo 

where M is the total mass initially. So that the region {(r, r) : < r < 6(t),t > 0} under 
consideration is transformed into the region {(x,t) : < x < M,t > 0}. 

Under the coordinate transformation (II. 191) . the equations (jl.6p - (|1.9p are transformed into 

d t p(x,t) = -p 2 d x (r n - 1 u), 

d t u(x,t) = r n ~ l {d x [p(* + 1p)dx(r n - l u) - P] - 2(n - l)fd x p} - f(x,r,t), (1.20) 
r n (x,t) = n Jq p~ 1 (y,t)dy, 

where (x,t) S (0, M) x (0,oo), with the initial data 

(p,u)\ t=Q = (p ,u ){x),r\ t=Q = r (x) = J q Po 1 (y)dy^j , (1-21) 
and the boundary conditions: 

u(0,t) = 0, (1.22) 
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{P - p(X + 2fj,)d x (r n - 1 u) + 2(n - %-} 



ic=M 



P r , t > 0. 



(1.23) 



It is standard that if we can solve the problem (jl.20p - (jl.23p . then the free boundary problem 
(fl~T|) - (fl~4"j) has a solution. 

From (|1.14p - (jl.l6p . it is easy to see that Poo(x) is the solution to the stationary system, 



Ar ™ X {pL)x = -/oo(x,roo), r%c(x) = n pj(y)dy, x G (0,M), 

io 



(1.24) 



(-U)- ('^ 



11.25) 



The results in [6J [18] show that the compressible Navier-Stokes system with the constant 
viscosity coefficient have the singularity at the vacuum. Considering the modified Navier-Stokes 
system in which the viscosity coefficient depends on the density, Liu, Xin and Yang in [9j proved 
that such system is local well-posedness. It is motivated by the physical consideration that in the 
derivation of the Navier-Stokes equations from the Boltzmann equation through the Chapman- 
Enskog expansion to the second order, cf. [1], the viscosity coefficient is a function of the 
temperature. If we consider the case of isentropic fluids, this dependence is reduced to the 
dependence on the density function. 

Since n > 2 and the viscosity coefficient p depends on p, the nonlinear term 2(n— l)^ud x p in 
(11.200 9 makes the analysis significantly different from the one-dimensional case [9"1 114 1 [TTl 119 1 120] . 
Considering the compressible spherically symmetric Navier-Stokes equations without a solid core, 
the techniques in the case of similar system with a solid core [U El [Til El EI] failed to be of use 
in our case, so we need obtain some new a priori estimates. 

For spherically symmetric solutions of the Navier-Stokes equations with constant viscosity, in 
[7], the author gave an information near the origin that the solution may develop vacuum region 
about the origin. The difficulty of this problem is to obtain the lower bound of the density p and 
the upper bound of the term -u. When the initial data are small in some sense, using some new 
a priori estimates on the solution, we can obtain the lower bound of the density and the upper 
bound of the term -u. The key ideas are using the classical continuity method and the result of 
Claim 1. In Claim 1, we want to prove that there is a small positive constant e\, such that, 
for any T > 0, if 

u 

I(t) = \\p(;t)- Poo \\ L ~+ -(;t) <2ei, Vt G [0,T], 

r l°° 



then 
Let 



I(t) < ei, V t G [0,T]. 



B[p,u,r] 



M 



(p-Poo) 2 +r 2n ~ 2+a (p-p 00 )l + 



+r 2n-2 n 2 + r 2n^+a {() l + B {r n-l u)x) , 



dx, 



where a = | — n. In Lemmas I3.3H3.81 we get some uniform a priori estimates (with respect to 
time) on the solution in the weighted Sobolev space and the upper bound of B[p,u,r]. Using 
the bound of B[p,u,r] and Sobolev's embedding Theorem, we can finish the proof of Claim 
1. Then, we will construct a weak solution by using the finite difference approximation. Our 
results show that: such system does not develop vacuum states or concentration states for all 
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time, and the interface dQ T propagates with finite speed. Since these estimates of the solution 
are uniform in time, we could show that the solution converges to a stationary one as time tends 
to infinity. Moreover, we construct various Lyapunov functionals and obtained the stabilization 
rate estimates of exponential type. 

We now briefly review the previous works in this direction. For the related free boundary 
problem of one-dimensional isentropic fluids with density-dependent viscosity (like p(p) = cp e ), 
see [HI \T7\ \19\ 120] and the references therein. For the spherically symmetric solutions of the 
Navier-Stokes equations with a free boundary, see [TJ [2l [Til Eil EI] etc.. Ducomet-Zlotnik|2| [2~T] 
studied the similar system with a solid core and without the nonlinear term 2(n — l)^ud x p. Also 
see Lions[8] and Vaigant-Kazhikhov|16j for multidimensional isentropic fluids. For the related 
stabilization rate estimates in the one-dimensional case, see [3| 110 } [121 1151 120] etc.. 

Main assumptions on ci, C2, and 7 can be stated as follows: 

(Al) condition ([TT7]) or ([Tig]) holds; 
(A2) 6 > 0. c\ and c 2 satisfy that 



Under the above assumptions (A1)-(A2), we will prove the existence of a global weak solution 
to the initial-boundary value problem (jl.2U|) - (jl,23|) in the sense of the following definition. 

Definition 1.1. A pair of functions (p,u,r)(x,t) is called a global weak solution to the initial 
boundary value problem (|1.20j) - (|1.23j) . if for any T > 0, 



d > 0, 2ci + nc 2 > 



and 



[2cjq + c 2 (2n - 2 + a)] 2 - 4(2ci + c 2 )[2 Cl (n - 1) + c 2 (n - l)(n - 1 + a)} < 0, (1.26) 
where a = | — n. 



p,u e L°°([0,M] x [0,T])nC 1 ([0,T];L 2 ([0,M])), 
r GC 1 ([0,T];L oo ([0,M])), 
(r n ~ 2 u) x , (r n ~ l ) x G L oo ([0,T];L ri_ ^([0,M])), 



and 



(r"-^), G L°°([0, M] x [0, T]) n CJ ([0, T]; L 2 ([0, M])). 



Furthermore, the following equations hold: 



Pi + P 2 (r n 1 u) x = 0, p(x, 0) = p (x) a.e. 




= n 




and 




/ / [mp t + (P-p(\ + 2p)(r n - 1 u) x )(r n - 1 ^) 



X 



+2(n - l)p{r n ~ 2 u^) x - f(x, r, t)tp]dxdt 




Pr{r n ~ l il))(M,t)dt - / u (x)ip(x, 0)dx, 



for any test function ip(x,t) G C^°(Q) with O = {(x,i) : < x < M, i > 0}. 
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In what follows, we always use C{Ci) to denote a generic positive constant depending only 
on the initial data, independent of the given time T. 

We now state the main theorems in this paper. Let p = min xg r 0i M] Poo and p = max xg j 0; Af] Poo- 
Theorem 1.1. Under the conditions M.lUty - kl.lty) and [Al]-[A2], there exists a positive constant 
eq > 0, such that if 

IK/!, AP)|| LoonL i + ||(AP)'|| i2 + ||/ 2 || i2 < eo, (1-27) 
IlPo - Poolll- + B[po,u ,r } < e§, (1.28) 
then the system il.20\) - [1.23\) has a unique global weak solution (p,u,r) satisfying 



p(x,t) G 



1 3_ 
2* 2~ P 



u , 



r n (x,t) G [C- X x,Cx], 
<C\\d x (r n -\){;t)\\L~ <Ce , 



B[p,u,r]<Ce 2 , 
for all t > and x G [0, M]. Furthermore, we have 

/o 



lim 



X n u + X~ 



(p% - (pL), 



dx = 0, 



lim ||u(-,t)|| L oo + \\p(-,t) - Poo(-)IU°° + \\r(;t) - r^-)]]^ 

t— >+oo 

Remark 1.1. In fact, assumption (|1.26p is give a restriction on -, i.e. 

-18 + 8n + 8n 2 - 8\/3(n - \)y/An - 3 
9 - 12n + 4n 2 
-18 + 8n + 8n 2 + 8>/3(ti - l)V4n-3 
9 - 12n + 4n 2 ' 



0. 



A 

< - < 
P 



(1.29) 

(1.30) 
(1.31) 



If n = 3, we can choose J = fa. g (|( 13 _ |( 13 + 8^3)). 

Remark 1.2. We can choose the constant eo as in (|3.75p . 

The proof of the uniqueness part of Theorem 11.11 also shows that the continuous dependence 
of the solution on the initial data holds. We may state the following result without a proof. 

Theorem 1.2. For each i = 1,2, let (pi,Ui,ri) be the solution to the system il.20\) - [1.23\) with 
the initial data (poi,uoi> r Oi)> which satisfy regularity conditions $1.29\) - [l.31\) . Then, we have 



M 



[(ui - u 2 ) 2 + (pi - p 2 ) 2 + x n(n- r 2 ) 2 ](x, t)dx 
-M 



< Ce 



ct 



[(«oi - u 02 ) + (poi - P02) + x » (r i - r 02 ) }dx, 







for all t > 0. 
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Theorem 1.3. Under the assumptions of Theorem \l.l\ and 



hit) + f 2 {t) + |AP(t)| + |(AP)'(t)| < Ce~ a »\ (1.32) 



where ao is a positive constant, then we have 



J {r 2n ' 2+a (p - Poo f x + r 2n - 2+a [d x (p l+e d x (r n - l u))) 2 + r a u 2 } dx < Ce~ 

(-,{r n - 1 u) x ) (;t) + \\p(.,t) -poc(-)IU- + IK,*) -rooOlk- < Ce" 
/or i > ; where a is a positive constant. 



at 



at 



Remark 1.3. Considering the general case that (/i,A)(-) G C(M+) H W, ' ( f°(]R+), under the as- 
sumptions (fITUD - (frT2"| ). (Al) and 

p{p) > 0, 2p(p) + n\(p) > 0, 

[2pa + A(2?i - 2 + a)] 2 - 4(2p + A)[2/x(n - 1) + X(n - l)(n - 1 + a)] < 0, 
for all p € [gP, §p], we can obtain the same results. 

Remark 1.4. In this paper, we study the case of 7 > 1 and prove the main results in this case 
only, since the case of 7 = 1 can be discussed through the similar process. The main different is 
that (|2.10p is replaced by 

S[V} = J ( AlnV x + PooVx + J Gx(nh)^dh) dx, 
when 7 = 1 and n = 2. 

The rest of this paper is organized as follows. First, we obtain the existence and uniqueness 
of the solution to the stationary problem in Section [2j In Section [3l we will prove some a priori 
estimates which will be used to obtain global existence of the weak solutions. In Section HI using 
the finite difference approximation and a priori estimates obtained in Section [31 we prove the 
existence part of Theorem 11.11 In Section \5\ we will prove the uniqueness of the weak solution. 
In Section [6l We show that the solution of the free boundary problem tends to a stationary one, 
as t — > +00. In Section [3 we will obtain the stabilization rate estimates of exponential type on 
the solution by constructing some Lyapunov functionals. 



2 The stationary problem 

We start with a proof of the existence of a positive solution to the Lagrangian stationary problem. 
Zlotnik-Ducomet [2 1 j studied the stationary problem with a solid core r > ro > 0. Using similar 
arguments as that in [21] , we can obtain the following results for the stationary problem without 
a solid core. 



Proposition 2.1. // 



or 



2n - 2 , . 

7> 2.1 

n 

2n — 2 , „ 2-n 2 

7 = and Gn—Mn < 2A, (2.2) 

n 
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or 

2n - 2 G 2- n 2 Ji^l 

0<7< — - — and P Xl + —n— M~5 z n < A5j, (2.3) 

n 

( 2n-2 \ 2n-2-n7 

— — — ^—3- ) , then the Lagrangian stationary problem U.2A\)- [1.25\) has a 

(n-l)GMn J 

positive solution p^ G W^QO, M]), where (3 G [1, rz^) is a constant. 
Proof. We introduce the nonlinear operator 

I : K —> W^([0,M]), 

where K = if £ C([0, Ml) : min fix) > ( — ^ ) 7 L by setting 

xe[0,M] V ^ / 



'(/)(*) 



( Poc + f x M G^^dy 



V 



,4 



with ry(x) = n Jq f 1 (y)dy, x G [0,M]. We can restate the problem (|1.24|) - (|1.25|) as the fixed- 
point problem 

Poo = I(poo)- (2.4) 



For all / G K s = {/ G K : f < 5} with 6 > (^) ' , we have 
and 

, , lv „ 2n — 2 2n — 2 / 2— n 

Poo<A{I{f)f < Pco + GS—n-— x~dx 

Jo 

G 2n — 2 2— ti 2_ 
= Poo + —5 n n n Mn. 

2n-2 2 — n 2 

If 7 > 2n n " 2 , then /(if^ ) c Kg x , where S\ is a positive constant satisfying P 00 +^S 1 n n~^~ Mn < 
A5j. And one can immediately verify that I is a compact operator on Ks 1 ■ Since K$ 1 is a convex 
closed bounded non-empty subset of C([0, M]), the problem (|2.4j) has a solution p G iQj by 
Schauder's fixed point theorem. 

If 7 = 2n ~ 2 and Gn~^~ M~ < 2 A, then I(K$ 2 ) C -K"<5 2 , where 82 is a positive constant 

In — 1 2 — n 2 

satisfying P oc + ^5 2 n n~ Mn < A5j. 
If 7 < ^ and 



2n-2 \ 2 n-2-ri7 / « 2n ~ 2 \ 2n-2-n-y 



Poo + — n~Mn [ £ I < A ' 



2 \ (n - 1)GM» / \ (n - l)GAf s / 

then I(Ks 3 ) C ii^, where 

\(n-l)GM» / 

We can finish the proof of the theorem immediately. □ 
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Letting Vx = using the equalit y = (V^x, one can eliminate the function poo from the 
Lagrangian stationary problem (|1.24p - (|1.25p and obtain an equivalent boundary- value problem 
for a non-linear second-order ODE: 

2-2n 2 ~ 2 " 

{A{V OQ )-^) x = -Gxn—Voc" , x G (0,M), (2.5) 

Foo(0) = 0, {VMM) = (j^-j 7 , (2.6) 

for a function V^ G C X {[Q,M]) such that (Voo) x > 0. 

In accordance with the method of small perturbations, we replace Voo by V = Voo + W with 
small W and linearized the operator in the last problem: 

, „/, 2-2n 2-2n 

{A{V)~ 1 ) X + Gxn— V— 
= (-T^VooJr -1 ^)* + ( 2 " 2n)Gx(nF 00 )^r i W + . . . , x G (0, M), 

f(o) = o + w(o), A(v x yy\ x=M -Poo = -iahv^-^w^m + ..., 

up to the terms of the second order of smallness with respect to W. We define the linearized 
operator 

L[W] = {- 1 Ap^ 1 W x )x + (2 - 2n)Gx(nV 00 ) 2 -^W, W G K Q , (2.7) 
where K = {W £ C x ([0, M]) : W(0) = 0, W X (M) = 0}. It is easy to get 

(L[W],W) = J (7^(Poo) 1+7 W^ - (2n - 2)Gx(nV r 00 ) 2: ^ !i iy 2 ^ ds, 17 G if . 

Let 

r-AI 

J[W]:=J (- f A(p 0C ) 1+ ^W^-(2n-2)Gx(nV 00 )^ 1 W 2 )dx, (2.8) 

for W G Jfi = {/ G ^([O, M]) : /(0) = 0}. 

We say a stationary solution is statically stable if 

> 5 3 (\\Wx(x)\\l H0M) + \\ X - l W(x)\\ 

1 2 (0,M)J > ( 2 - 9 ) 

for some £3 > and all W £ K\. 

Now, the static potential energy takes the following form: 

S[V] = J (^3l^ 1 " 7 + Po °^ + i / 1 Gx(nh)^dh\dx. (2.10) 

We call 7 G = {/ G C 1 ([0, M]) : /(0) = 0, min(/ x ) > 0} is a point of loca/ quadratic minimum 
ofS if 

5 [V + W ] _ > «j 4 (\\W x (x)f L2{0M) + Wx^WWW^qm) > (2-H) 

for all W £ K\ and ||W||cfi([o,Af]) — <^5> f° r some ^ > and <5s > 0. 

We can clarify the variational sense of the definition of statically stable as follows. 

Proposition 2.2. A function V G K% is a point of local quadratic minimum of S if and only if 
V = Voo is a solution of the problem \2. 5\) - [2~6)) and satisfies static stability condition A2.9\) . 



9 



Proof. Let V £ K2, W £ K\ and || W||ci([o,A/]) = 1- Using Taylor's formula, we have 
S [V + eW] = S[V] + SS[V](eW) + \^S[V + reW]\ T=f , 

where 

f M / 2-2n \ 

6S[V](eW) = J [-A{V x )^eW x + P 00 tW x + Gx(nV) — eWj dx, 



and 

dr 2 



S[V + reW] = / (-fA^ + reW^-^eW,) 2 
Jo 

-(2n - 2)Gx(n(y + reW))^ {eW) 2 ) 
for all |e| < mi ^ v and some f € [0,1]. If (|2.1ip holds, we have 

' /2 5[y + T6iy] < ce 2 (||w x (x)||| 2(0 » /} + \\ x - l w{ x )\\l H 



Ticrr|Nw<; ^l"x^;|| L 2 (0iM ) -r 11^ kk V*-;ilZ,2(0,M) 

and 

Ce 2 (||^(x)||i 2(0>M) + II^WCx)!^^) + e<5S[y](W0 > 
for all |e| G (0,min((55, mi „y )) and llWllcir^M]) = !• Thus, we obtain 

5S[V]{W) = 0, 



i.e. 



J {-A{V X )^W X + Poo^ + Gx(nV)^V) dx = 0, 



for all W £ i^x and ||W||ciro m] = 1> that is, V is a stationary point of S and a solution of the 

d 2 

dr- 

d 2 



problem Q2.5| ) - (|2.6p . We can rewrite 4^S[V + reW] as follows 



, 9 -S[V + reW] = 5 2 S[V]{eW) + Si, 
dr z 



where 5 2 S[V]{eW) = £iS[V + reW]\ T=0 and 



d 2 

—^S[V + reW] - 5 2 S[V](eW) 
dr z 



< Ce (\\eW x (x)\\ 2 L2[[0M]) + \\x- l eW(x)\\ 2 LH[om ) . 



Thus, we obtain 

5 2 S[V](eW) > (<5 4 - Ce) (\\eW x (x)\\ 2 LH[0M]) + Wx^eW^W 

l 2 ([0,Af]) 

for all e G (0,min(<5 5 , 5^4, J^)) and ||W||c :l ([o,M]) = 1- Moreover, we have 

J[W] := 5 2 S[V](W) > | (||^(x)||| 2([0iM]) + llx-M^lllw,])) , (2-12) 
for all W E K\. 

If V = Vrx, is a solution of the problem (|2.5p - (|2.6p and satisfies static stability condition (|2.9p . 
we can prove Voo is a point of local quadratic minimum of P easily. □ 
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Proposition 2.3. IfV = V QO is a solution of the problem \2.5\) -\2. 6 \) and 7 > 2 ^- 2 -, then \2. 9\) 
and 12~T1\) hold. 

2 — 2n 

Proof. From (Ap2o) x = —G = — Gx(nVoo)~^~ , using integration by parts, we have 
= J ( Kl A{ Poo ) l+ ^W%-{2n-2)Gx{nV 00 ) 2 -^W 2 ^dx 



M 



(7^(fe) 1+1 ^ 2 + (2n - 2)^(/>2 o ) ;E (ny 0O )^ 1 ^ 2 ) cfe 



/ (7^(Poa) 1+7 W* - 2(2n - 2)^(11^) - 1 WW, 
io 

+ Ap£ l V~ 2 W 2 )dx + (2n - 2)P 00 (M) 

/ W 2 \ 

J W + 2n - 2)Poo — — (M), for all W e K - (2.13) 



If 7 > 2n^2 h 



M 2n-2 . ^ / W x2 



W > / ^L +7 W x — dx. (2.14) 

Jo n V P00VW 

If (|2.9p not holds, we have for any integer m > 1, there exists W m £ Kq and HW^IIcirm.jWl) = 1 
such that 

J[W m ] < ^ (\\(W m Ux)f L2{0>M) + \\x- l W m {x)f L ^ M) ) . (2.15) 
Then, there is a subsequence m —* 00 for which 

W m ^I¥ in C([0,M]), 
(W / m ) x — 51 Wj; in L 2 ([0,M]). 

From ([2TT3>([2TT5]) . we have 

w x = -^ r , xe(o,M), 

Poo 'CO 

and W(0) = W(M) = 0. Thus, we obtain W = 0. It is a contradiction. 

Therefore, if 7 > 2 ^, then (|2U|) holds. From Proposition we can obtain (|2TTT|) 

immediately. □ 

Now, we shall use the shooting method to prove the uniqueness of the solution. 

Proposition 2.4. Under the assumptions \2.1\) - l2~3ty , the Lagrangian stationary problem [1.24 )- 
il.25]) has a unique positive solution p^. 

Proof. We consider the Cauchy problem 

(Ap2o) x = -GxinV^) 2 ^ , (V^ = p£, x e (0,M), (2.16) 
P^\ x=0 = °, ^L =0 = 0, (2.17) 
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for the unknown functions poo(a,x) and V oa (a,x), where a > is the shooting parameter. For 
each a > 0, using similar arguments as that in Proposition 12.11 we can obtain the existence of 
the solution to this problem, satisfying 

(j oc 2x 

Poo((T,x) G [~,ff], Voofax) G [-,—], xG [0,M„], (2.18) 

2 GO 

Poo G iy 1,/3 ([0, Mo]), G C\[0, M ]), (2.19) 

2n — 2 2 — ra — , \ 

where Mo is a positive constant satisfying Aa' — a n ^-n n M ™ > A (?) and Mo < M. If 
there exist two solutions (pi, V\) and (p2, V2) to this problem satisfying 

Pi €W 1 ' /3 ([0,M < ]), xG [0,M<], (2.20) 

where Mi G (0, M], i = 1, 2. From (I2.20|) . there exists a positive constant M 3 G (0, min{Mi, M 2 }) 
such that 

(j x 2x 

Pi(x) G [-, a] and Vi{x) G [-, — ], x G [0, M 3 ], i = 1,2. 

2 a a 

Then, we have 

C X 2-2n 2 - 2n 2-2n [X 2 _ 2n [V 

Ap\-Ap\ = \ Gyn—(V 2 » -V x n )dy<C y~ / K - P; r 



and 



2 



IIPI - H|L°°([0,e]) < C||Pl - HU°°([0,£]) / ?/ n ^ < C<t£™ - P2||L°°([0,el), 

Jo 

for all x, e G (0, M3]. Choosing e < Co- 2 , we have 

Pi = P2, for all x G [0, e]. 

Considering the Cauchy problem 

{A P 1) X = -GxinVoof-^ 1 , (Voo)x = p^, 1 , x G (|,Af), (2.21) 

Poo| x . = | = Pi(cr, -), Kxd| x= | = y Pi {cr,y)dy, (2.22) 

using the classical ODE theory, we have pi(x) = p2(x), x G [|, min{Mi, M2}]. Thus, for each 
a > 0, there exists a unique solution to the problem (|2.16p - (|2.17p satisfying poo(x, a) > for 
x G [0, M CT ), where either Poo\ x=M = and M CT G (0, M) or M CT = M. 

Clearly, if p^ is a solution to the problem (jl.24j) - (jl.25j) . then p^ satisfying (I2.16p - (l2.17p for 
some a > 0. We will show that this can be possible only for one value of a. Using similar 
arguments as that in the above part and in [5] (§V.3), we obtain that (d a p2o, d^Voo) is well 
defined and satisfies the linear Cauchy problem 

Md,plo)x = (2n - 2)Gx(nV 00 )^ 1 d (T V 00 , (d^V^ = --p" 7 " 1 ^, (2.23) 

7 

where x G [0, M a ), 

d CT p£4 =0 = l, 0^1^ = 0. (2.24) 

It is easy to see that 

9 ct p2o > 0, {dvV^ < 0, dcVoo < 
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hold on [0,M 4 ), where either d a p2o\ x=M4 = and M 4 G (0,M CT ) or M 4 = M CT . We claim that 
only M 4 = M CT can occur. 

Assume that M 4 G (0,Af<r). Letting cj> = Ap^d^V^x + ^2 A 9 a p2o(Voo) x , from (j2"T6D and 
j223D, we have 



M 4 







In — A 



AU 



By the estimates pooi&jM^) > 0, d a p2o\ x=M = 0, <9 cr V r 0O | a;=M4 < and the initial condition 
(f2T7l) and OTTHh . we get 



M 4 

(Adx < 0. 





On the other hand, from ()2.16|) and (12. 23ft . we have 

4> = Ap" 1 ^^-^— ~ -) > 0, x G (0,M 4 ). 
in — l 7 

It is a contradiction. 
Thus, we obtain 

Poo > 0, d a poo > 0, x G (0,Af CT ). 

and M CT is non-decreasing on c G (0, 00). Therefore, for each fixed point x G [0, sup cr>0 M a ), 

1 

the function poo(<7, 2;) is strictly increasing on a > (^j 2 -) 7 , and satisfying Ap2o\ _ M = Poo for at 
most one value of a. □ 

Using the properties of the transformation (|1.19p and Propositions 12. 1112.41 we can obtain the 
following proposition immediately. 

Proposition 2.5. Under the assumptions jl2.1\ )- (2Jfy . the Eulerian stationary problem |777^p- 
il.16]) has a unique positive solution (poo> loo), satisfying < p < Poo(r) < p < 00, (poo)r( r ) < 0, 
< r < Zqo with loo < +00. 

3 A priori estimates 

From (jl.lip . (|1.20p and (|1.24j) . we could obtain the following lemma easily. 
Lemma 3.1. Under the assumptions of Theorem \l.l\ we have 

n = u, (3.1) 
rM f< 

ApUx) = Poo+ ^d—dy, (3.2) 
Jx Too (y) 

^) 7 < Poo < p < 00, r» (x) G [C-'x, Cx], (3-3) 



d , . „ , ^ „ x 
dx 



(ApUx)) = -G-^, (3.4) 



for allxe [0,M]. 
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Lemma 3.2. Under the assumptions of Theorem \l.l[ we have 



dtj \2 7-1 p J 1 a»-i J 

+ /"{(-el + csJ^Kr^u),] 2 + ^-^-c lP ^(r n ^u x - -f}dx 
Jq n n rp 



/ Afudx - Ai>(ur n-1 )(M,t). 
•/ o 



(3.5) 



Proof. Multiplying (|1.2(jp 9 by u, integrating the resulting equation over [0, M], using integration 
by parts and the boundary conditions (|1.22j) - f|1.23p . we obtain 



— / -u 2 dx- / Ap^d x {r n - l u)dx 
dt Jo 2 Jq 

-M 



+ J {(2 Cl + c 2 )p 1+e [(r n - l u) x ] 2 - 2d(n - l)p e {r n ~ 2 u 2 ) x ) dx 



-P T {ur n -"){M,t)- / fudx. 

'o 



From (jl.20p . we have 



/ Ap<d x {r n - l u)dx = ^ / rp^dx, 

Jo dt Jo 7 - 1 



(3.6) 



(3.7) 



and 



-P r (ur n ^)(M,t) = -P 00 (r 4 r n " 1 )(M,t) - AP(ur n " 1 )(M, t 

-M 



d r M p 

— / — dx- AP(nr"- 1 )(M,i), 



fudx 



d f M f T x f M 
dtl h G ^ dsdX ~J Afudx > 



(2 Cl + c 2 )p l+e (r n - l u)l - 2d(n - l)p b '(r"~V) 



= {-c l+ c 2 )p^{r^uY x + 
n 

From (|3.6p - (|3.10p . we obtain (|3.5p immediately. 



1 l+0/'_n-l. 



u 



cip (r u x j 

rp 



(3.8) 
(3.9) 



(3.10) 
□ 



Claim 1: Under the assumptions of Theorem ll.il there is a small positive constant e±, such 
that, for any T > 0, if 



J(t) = \\ P (-,t) -pooh* + -(■,*) <2ei, Vi G [0,T], 



(3.11) 



then 



7(t) < d, Vi e [0,T]. 

Using the results in Lemmas I3.3ti3.8t we can give the definition of e\ in (|3.74p and finish the 
proof of Claim 1. 
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Lemma 3.3. Under the assumptions of Theorem \l.l\ and h3.11\) , if e\ is small enough, we obtain 



p(x,t) G 



1 3_ 
2* 2 P 



r n (x,t) G [C^x^Cx], 

u{;t)\\i* + \\p(;t) - Pooh* + \\r- n (r n - r^U* < do, 
t rM / „,2" 



(3.12) 

(3.13) 
(3.14) 

(3.15) 



j J (u 2 + r 2n - 2 u 2 x + dxds < a4, 
for all t G [0, T) and x G [0, M). 

Proof. From Lemma 13.11 and (|3.1ip . we can easily obtain the estimate (|3.12p when 2e\ < \p. 
From CEO} 3 and (|3TT2]) . we can obtain (|3TT3j) immediately. From ([2TTP]) . ([3T2"j) and (|3T5) . we have 



M 



/o I n 

-M 



l+eun-l^2 , 2(n 1) 1+e ^ _ u_ dx 
n rp 1 



/ A/udx - AP(w n_1 )(M, i) 

JO 



(3.16) 



where Voo = ^ and V = £. From (OTP . ([2TTH . (l3TT2|) - (l3TT3l) and Proposition EES we have 



C 



-i 



(p-Poo) + 



2 , 



21 



V 2 
' oo 



(P-Poo) + 



2 Ay- Voo? 



V 2 

v oo 



< 5[y] - sfvy < c 

JO 

when || V - V^Hci^m]) < C^l < *5> and 

|AP(nr n - 1 )(M,t)| < Ce ^(r^u^dx 
From (fL27j) - (fL28]) . <|3TT2|> and (|33B]) - ([3Tl8j) . we obtain 



dx, 



(3.17) 



(3.18) 



r-M 

/ {u 2 + (p- Poo? + r" 
Jo 



+ 



t r M 



J 



(r™"^) 2 + (r"- 1 ^ - — ) 2 \ dxds 



rp 



< Ce 2 + C fi(s)\\u(;s)\\ L 2ds, 



using Gronwall's inequality and (jl.27|) . we can obtain (|3.14n -f f3.15[) immediately. 
Lemma 3.4. Under the assumptions of Lemma \3.3\ if €q is small enough, we obtain 



t r M 



J 



[(p - Poo? + r" 2n (r" - r^) 2 ] dxds < C 3 e 2 , 



(3.19) 
□ 

(3.20) 



for allte [0,T\. 
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Proof. Multiplying (jl.20|) '? by r "(^ — integrating the resulting equation over [0, M], using 
integration by parts and the boundary conditions (|1.22p - (|1.23p . we obtain 



M 



n n 



M 



oo 



Ut ,r 



o 



r n 1 n 



v 



-)dx + AP 





' oo 1 




\ n 


n J 


x=M 







M 



M l-n 

A/ (r n - r^)dx + I 2d(n - l)p & ( -(- - ^ 



o 



r n 



v 



dx 



+ 



M 



(2 C1 + PaJp^flSrCr"" 1 ^^ 1 - P^dx 



i=l 



We can rewrite the left hand side of (|3.2ip as follows 



L.H.S of (I3T2TD 

M 



-1\2 



(2n - 2)Gxr 



2-3n / r 



+ 



91 (ft 



-i 



1 n n 



re n 



dx 



where 



and 



\9i\ 



A(p 



P 1 - Poo 1 



7^p 



1+7 

oo 



<C 4 ei 



152 1 



+ (2n - 2)Gxr 

n n 



2-n 
oo 



< C 4 €l. 



From (|2.9p . we have 

L.H.S of (13311 > (2C 5 - C 4 ei 
> C 5 



M 



M 



(p-'-p-Jf + r^ 

n n 



n 

„n \ 2' 



' oo \2 

n 



dx 



dx, 



(3.21) 



(3.22) 



when 6461 < C5. 

From (|3.11|) and (|3.12|) - (|3.13|) . using integration by parts, we can estimate Ii as follows. 

= 4 f » f£ - + l M j ( 1 + '" - 1 >- N ) d. 

dt Jq r n 1 \ n n J J \n nr n J 

< "4 / M -^(--^dx + C [ M u 2 dx, (3.23) 



Jq r n 1 \ n n 



AP J (p- 1 - p-J)dx <^J (P" 1 " P-Jfdx + C\AP\\ (3.24) 
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C l' M 



2n 



n n 



dx + Cfl 



(3.25) 



I A < 



n rM 
10 



(P^-Poo 1 )^ 



-2n 



n n 



dx 



+ °r { [{rn ~ lu)x]2+ r^) dX 



and 



h<-^J Q (P- 1 ~ P'Jfdx + C J (r^u&dx. 
From (^2T]) - (pT37) . we get 



(3.26) 
(3.27) 



d 



M 



U 



dt ./ nr 

-M 



M 



{r n _ r n )dx + C / [( -1 _ -1)2 + r -2n ( _ ^o)^ 

1 Jo n n 

2 s 



< ^ r 2n-2 u 2 + ^ da; + c ,(| A p|2 + / 2) 



And from (1332]) - (13351) . we obtain (I3T201) immediately. 



(3.28) 



□ 



From now on, we study the case of 8 > and prove the main results in this case only, since 
the case of 9 = can be discussed through the similar process. 



Lemma 3.5. Under the assumptions of Lemma \3.3[ if e\ is small enough, we obtain 



M rt r M 

2n-2/„ „ \2i/ .\ t . / / rJn-2/ 



[r Zn -'(p- Poo y x }(x,t)dx+ / [r Zn - 2 ( P - Poo y x ](x,s)dxds<C 6 4, 
o Jo Jo 



(3.29) 



for all t £ [0,T]. 

Proof. From (|1.20p . we have 



d t H + H 

2 Cl + c 2 

^^ p i-e u+ ( 2ci+c2 _ )(n _ ^ _ ( } 

2ci + c 2 6* 



(3.30) 



where -ff = n+ 2ci + C2 r" l [p e — / 0^ o ) x . Multiplying (|3.30|) by i?, integrating the resulting equation 
over [0, M], using the Cauchy-Schwarz inequality, we obtain 

d_ ? M 
dt 



< 



\ H 2 (x,t)dx + C 7 H 2 (x,t)dx 
Jo Jo 

C J {\Hp^- e u\ + \^H 2 \ + | ^-H\ + \AfH\j dx 



+ 



M 



,7— fl^ra— 1 



M 

\H\dx + C I {r^uip^xHldx 
o 
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< c 



M / n 4 x 2 r 2n-4 \ j /-AT 

u2 + ^2 + r 4n-4 u 2 J dx + {- + C 8 ei)C 7 I H 2 dx 







+C 



a/ 



G— -r + P „ a {Apl 



„n— 1 



-y-9 
Pic 



dx + Cfl 



From (pTij) and (f3TT2]) . we have 



Ai 



r n— 1 



7-0 
Pi) 



d.T 



A/ 



/•M 

< C [(r - r^) 2 + (p - p^) 2 ] dx. 

JO 

Then, if e x < 1 and C 8 ei < ±, from (f3T3ll - ([332j) . we obtain 

d f' M C f M 

— / H 2 (x,t)dx + -1 / H 2 (x,t)dx 

/■Af 

< C {u 2 + (r- roD ) 2 + (p- poo) 2 ) dx + C/ 2 . 
•/ 

From d5121) - d3151) . (13T20D and (ET551) . we obtain ([3T29D immediately. 

Lemma 3.6. Under the assumptions of Lemma \3.3l if e\ is small enough, we obtain 



X Q^^P 1 1 



„n— 1 



7-^„2n-2 
poo '00 



dx 



(p(M,t)- Poo (M)Y+ / (p(M, S )-p oo (M)) 2 d S <C 10 eg, 

Jo 



t r M 




JO 



r2 m (p- Poo ) 2 dxds < Cueg, 



Ai 



(ri-"V)(x,t)dx + T l M r\~ m (r 2n - 2 u 2 x + ^) dxds < C 12 e 2 , 

(r 2n-2 +5 - m(/) _ poo) 2 )dx + / / {r 2n-2+- 2 - m{p _ pgo) 2 )dxda < ^ 

11 io io 

/•Af 

\\p(-,t) - Poo(-)\\l°° + / \(P ~ Poo)x\(x,t)dx < Ci 6 e , 



kOM) - »"oo(«)| < Ci 6 e x«, x£[0,M], 
for all t G [0, T] and m = 0, 1, . . . ,n — 1. 

Proof. From (|1.2U|) i and the boundary condition (|1.23|) . we have 

(2ci + c 2 ) 



W)(M,t) = -2d(n- 1) (/^) (M,t). 



Multiplying the above equality by p 6 (M,t) — p 6 OQ {M), we obtain 
2ci + c 2 d 



26* di 



+ (p°(M,t) - p^(M)){Api(M,t) - P a 



(3.31) 



(3.32) 



(3.33) 
□ 

(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 
(3.39) 



x=M 



IS 



= -2d(n - 1) \ P e -(p e - p 6 J\ + AP(p d (M,t) - pl(M)). 
Combining (|3.12p - (|3.13p . using the Cauchy-Schwarz inequality, we get 

j t (p 9 (M, t) - pl{M)f + C-\p e (M, t) - pl(M)) 2 

r-M 

< C\AP\ 2 + C(u 2 r n )(M,t) =C\AP\ 2 + C d x (u 2 r n )dx 

Jo 

< C \AP\ 2 + C J\r 2n - 2 u 2 x + ^)dx. 



(3.40) 



Integrating the above inequality over [0, f], using the estimates (|3. 12)) and (j3. 15|) . we can obtain 
([3341) . 

From (pTTijl - ipOH] ). (pTZOj) and (^29|) . we know that the estimates (l3T35l) - (l3T37]) hold with 
m = 0. 

Claim 2: If that (l3T35l) - (l3T37l ) hold with m < k, k G [0, n-2], then the estimates (l3T35]) - (l3T37| ) 
hold with m = k + 1. 

We could prove Claim 2 as follows. Let a>k = \ — k — 1. Using Holder's inequality, we have 

-M 

r ak (p~ poofdx 

M / r-M \ 2 

r a * (p(M,s)- Poo (M)- d x (p- Poo )dy) dx 



o 



< C(p(M,s)- Poo{M)f 

r-M r-M r-M 

+C r ak r 2n ~ 2+ak -^{p- Poo )ldy r 2 - 2n - a ^dydx 



JO Jx 

r-M 

< C(p(M, s) - Poo(M)) 2 + C / r 2n ~ 2+a ^ (p - Poo ) 2 x dx. (3.41) 

Jo 

From d333J, ([337D(m = fc) and ([531]) , we obtain ([335]) (m = k + 1). 

Multiplying (|1.2U|) 9 by ur ak , integrating the resulting equation over [0, M], using integration 
by parts and the boundary conditions (|1.22j) - f|1.23j) . we obtain 

J r-M -i r-M 

4- / -r^u 2 dx- / ^r^-Wdx 
dtj 2 J 2 

M 

[(2ci + c 2 )p 1+e (r n - 1 u) x (r n ~ 1+ak u) x - 2 Cl (n - l)p e (r n - 2+ak u 2 ) x ]dx 

rM r-M 

+ / A{p"t - pl ){r n - 1+at >u) x dx - / Afur ak dx 



r-M r-M 

-AP(ur n - 1+ak )(M,t) + / (r n ~ 1+Q *u) x 

Jo Jx 



Gy Gy 



J2n-2 r 2n-2 
I oo ' 



dydx 



= J>i- (3.42) 
i=l 

We can estimate L\ as follows. 

-M 



-U = J {(2ci + c 2 ) 



r 2n-2 +akf) l + e u 2 
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+ [2a kCl + c 2 (2n - 2 + a k )]p d r n ~ 2+ak uu x 

+[2cj(n - 1) + c 2 (n - l)(n - 1 + afe)]p e-1 r a * _ V} da;. (3.43) 

Since 

[2acj + c 2 (2n - 2 + a)] 2 - 4(2cj + c 2 )[2ci(n - 1) + c 2 (n - l)(n - 1 + a)] < 0, 
where a = | — n and 

[c 2 (2n - 2)] 2 - 4(2ci + c 2 )[2ci(n - 1) + c 2 (n - l) 2 ] < 0, 

we have 

[2a k a + c 2 (2n - 2 + a fc )] 2 - 4(2ci + c 2 )[2c l {n - 1) + c 2 (n - l)(n - 1 + a k )] < 0. 
Then there exists a positive constant C13 such that 

- Li > C13 jf ^ (r 2 "- 2+Q V 1+e ^ + AV*-V) dx. (3.44) 
From (|3.12j) - (|3.13p . using the Cauchy-Schwarz inequality, we obtain 

"M rM 7 _ 7 



/•AZ /-M 7 _ 7 

= / A(p 7 - pl^r^+^uxdx + A(n-1 + a k )r a ^ l u^- ^dx 

Jo Jo P 

< J (r^-^p^ul + p^r^u^dx + C J r a * (p - Poo ) 2 dx, (3.45) 

1 r M 

L 3 < -C 13 / p^V^Vdx + Cfl (3.46) 
° Jo 



WW 

L 4 = - / AP<9 x (r"- 1+a ^)da; 

JO 

< ic 13 [ M ( r 2n - 2+a *p 1+9 ul + p^V^u 2 ) dx + C|AP| 2 (3.47) 







and 

*M 



L 5 < ic 13 J (r 2n - 2+ ^p 1+e u 2 x + /-V^- 2 n 2 ) dx 



< 



+C I /•'"■ 
Cl3 



M rM r Gy Gy , 

dy 



2n-2 ^2n-2 
x \ ' 00 



d.r 



(r^-2+a kp i+e u 2 +p e-l r a k -2 u 2^ dx + c r ak (p - Poo ) 2 dx. (3.48) 



' JO 

From flSZEED , (!3TT2l) - (l3TT3l and d33Z|> - ([53ED , we obtain 



-M r M / „, 2 - 

-•"- 2 n 2 + ^ ] d.r 



d r M r M / 

^- / (r a 'V)(M)dx + 2C7 14 1 / r afc r 2 
dt Jo Jo V 

< C(/ 2 + |AP| 2 ) + Ci 4 ei / r^— dx + C/ r a * (p - Poo ) 2 dx. 

Jo r Jo 
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When C± 4 ei < \, using the estimate (|3.35p (m = k + 1), we can get 



J fM r 

-I (r^u 2 )(x,t)dx + C^ J 

fM 

< C(fi+\AP\ 2 ) + C r a *{p- Poo fdx 

Jo 



a k / J2n-2 q 2 



ur 



r zn u z x + — ) dx 



(3.49) 



and (13361) (m = k + 1) holds. 
From (11.20p . we have 



9 t Fi + — ^—pi-OHx 
2ci + c 2 



2 2ci + c 2 P 



+ 



2ci + C2 a fc 

# (n - 1 + y) - 2ci(n - 1) 



(n-l + a)( 2ci + C2 ) n _ 2+ ^ 



"OoJs - f(x,r,t)r 2 



(3.50) 



where H\ = r £ u+ 2ci + C2 r n 1+ 2 fc (p e — p^x- Multiplying (j3.50p by J3i, integrating the resulting 
equation over [0, M], using the Cauchy-Schwarz inequality, we obtain 



d f M f M 
— / Hf(x,t)dx + C ls / Hf(x,t)dx 
dt Jo Jo 



M 







< C I [^H^-Vr— u\ + 



Hf + r—^u 2 ^ + r~?\AfHi\ ) dx 



+C 



< C 



+ r 2 



Hr rfl— 6 r n— 1 

+ ^-(^3c). 



-n— 1 



7" 
Poo 



|#i|cfe 



2 2n-4+a fc 



+(^ + Ci 7 ei)Ci 8 / H 2 dx + C 



M 



4n-4 
' oo 



u 2 )dx + C/f 



,7— 9 r n— 1 



„n— 1 



+ 



7—6 
Poo 



iP 



dx. (3.51) 



If e\ < l,Cne\ < 4, using the estimates (|3.1ip and ()3.12p - (|3.13p . we have 



J fM (~i fM 

-J H 2 (x,t)dx + ^J H 2 (x,t)dx 

pM 

< C / [r ak u 2 + r a " (p - Poo ) 2 + r a " (r - r^) 2 } dx + Cf 

Jo 

,-M 

< C [r ak u 2 + r a *{p- Poo ) 2 ]dx + Cfl 

Jo 



(3.52) 



And from (13.35p - (13.36p (m = k + 1), we have 



M ft fM 

Hf(x,t)dx+ / H 2 (x,s)dxds < Ce, 
Jo Jo 



From (I37T21) and (l3T35|) - ([3T36ll (m = k + 1), we obtain (13371) (m = k + 1) immediately and 
finish the proof of Claim 2. 
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From Claim 2, we obtain that the estimates (|3.35j) - f)3.37[) (m = 0, . . . , n — 1) hold. From 
([3.13P and (|3.37p . using Holder's inequality, we obtain 



M 



\(p-Pao)x\dx<U r 2n - 2+a (p- Poo y x dx 



M 



,-2n+2-a r 



dx < Ce . (3.53) 



From (|3.13p - (|3.14p and (|3.53p . using Sobolev's embedding Theorem, we could obtain (|3.38j) - (|3.39p 
immediately. □ 

Lemma 3.7. Under the assumptions of Lemma \3.S\ if ei is small enough, we obtain 

[ M (jt + r2n ~ 2u l) (x,t)dx + uUx,s)dxds<C 9 e 2 Q (l + \\(r n - 1 u) x \\ Lr J, (3.54) 

for all t G [0,T]. 



Proof. Multiplying (|1.2U|) 9 by Ut, integrating the resulting equation over [0, M], using integration 
by parts and the boundary conditions (|1.22p - (|1.23p . we obtain 



M pM 

u 2 dx + / (2c x + c 2 )p 1+9 (r n - 1 u) x (r n - 1 u t ) x dx 
JO 

M 

Ap<{r n ~ x u t ) x dx - P v {r n - x ut){M,t) 
o 

r-M rM 

+ / 2c\{n — l)p e (r n ~ 2 uut) x dx — / futdx 
Jo Jo 

4 



8=1 



From (|3.1ip . (|3.12p and (|3.15p . using the Cauchy-Schwarz inequality, we obtain 

r-M 

(2c l +c 2 )p 1+e (r n ' 1 u) x (r n - 1 u t ) x dx 
M 







dt ./o 



(2d + c 2 )(n - l)p l+e {r n - l u) x {r n - 2 u 2 ) x dx 



o 



!- / M (2C1 +C2) (1 + e)p 2 + d [{r^u), x fdx 
lo 1 

d_ [ M 2£l±^2 p l +d[{r n-l u)x] 2 dx 



dt 



o 



C(\\(r n -'u) x \\l 2 + J )(l + ||( r «-i n)x ||^J, 



LI' 



(3.55) 



(3.56) 



J rM rM 

- Ap'?(r n - 1 u) x dx+ A 1 p~< +1 [{r n - 1 u) x ] 2 dx 
dt Jo Jo 



M 



2A(n - l)p^-(r n - l u) 
o r 



f M _ u 2 

x dx + / Anin — l)/? 7 " 1 — ~dx 
Jo ' r 2 
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d f M 11 2 

< ± Ap\r^u) x dx + C{\\{r^u) x \\h+ ~ ! 



(3.57) 



iV 2 



d_ 

dt 



Jo 

i-M pM 

+ / (n - ^(Poo + AP)(r n - 2 u 2 ) x dxds + (AP)' / (r^^^yds 
Jo Jo 

d r M 



+ 



Li 



+ I(AP) 



/|2 



(3.58) 



and 



- / ci(n - lJp^r^-Vj.dx + / 20d(™ - lJ^-Kr"- 1 ^]^ 
at Jo Jo r 

fM u 2 t-M u 3 

- 8c 1 n(n-l)p 9 — (r n - 1 u) x dx+ 2n Cl (n - l)(n - 2)p 6 ~ l ^dx 
Jo r Jo r 



M 







3ci(n - l)(n - 2)p t) ^(r n - 1 u) x dx 



d 



M 



< - j o c l (n-l)p t) {r n - z u%dx + C{\\{r n - l u) x \\i l + 



Li 



N A <-^- / G^-dx + / (1 - n)Gxr- n u 2 dx + - / u 2 dx + Cj\. 
dt J r n 1 J 2 J 

From (^551) - (1330]) . using the fact that 



"M 



"M 



M 



M 



-(2 Cl + d)^^ 1 ^] 2 - d(n - l)p e {r n - 2 u 2 ) x dx 



1,2 



(n-1) 



2 v n 



7? 



-(-n +c 2 )p 1+9 [(r"- 1 n) :c ] 2 + i^—^c!/^- 1 ^ - ^) 2 }> ,/,-. 



we have 



M 



-ui dx 



+ 



d f M f 1 , 2 



eft 



2 v n 



u 



< 



d r r M r /* M 



'-c lP \r n -'u x - -Y dx 

n rp J 

Gy Gy 



Jln-2 r 2n-2 



dy 



(3.59) 



(3.60) 



dx 



+C{l + \\{r n -'u) x \\ L ~)(\\(r n -'u) x \\l 2 + 



Li 



+ /i 2 + l(AP) 



/\2\ 



(3.61) 



Integrating (|3.6ip over [0, t], using the estimates (|3.13p - (|3.15p and the Cauchy-Schwarz inequality, 
we can obtain (|3.54p . □ 



23 



Lemma 3.8. Under the assumptions of Lemma \3.3l we obtain 



M 



t rM 



r a u 2 )(x,t)dx + 




r 



2n-2+a 2 , «-2 2 



u xt + 7,0 u t) dxds < Cige , 



J 



■(■,t) +\\(r n - 1 uU;t)\\ LO o <C 20 e , 

f M fu 2 \ r* r M 

J l^2+r 2n - 2 u 2 x j (x,t)dx + J J u 2 (x,s)dxds <C 21 e 2 , 



(3.62) 
(3.63) 
(3.64) 



for all t G [0, T], where a 



n. 



Proof. We differentiate the equation (|1.20j) 9 with respect to t, multiply it by utr a and integrate 
it over [0,M], using the boundary conditions (ll,22|) - (ll.23p . then derive 



d_ 

dt 



1 



r a u 2 dx 



a 



M 



r a l uu 2 dx 



M r 



(2 C1 + c 2 ) P 1+e (r n - 1 u) x - Ap~< + Poo - 2ci(n - l)p l 



u 



x((n- l)r n - 2+a uu t ) x dx 



M 



f 







(2c l +c 2 )p l+e (r n ~ 1 u) x -Ap^ 



+ApZ -2c l {n-l)p°- (r n -' +a u t ) x dx 
r J 

rM 



cM 



f t r a u t dx 



+ / 2 Cl {n-l)d t {r n - l p e d x {-))r a u t dxds- 
Jo r Jo 

- [(n - l)AP(r n ~ 2+a uu t )(M,t) + (AP)' (r n - 1+a u t )(M,t)] 
:= J! + J 2 + J 3 + J 4 + J 5 - 

From (|3.12p . (|3.35f) and (I3.36|) . using the Cauchy-Schwarz inequality, we obtain 



(3.65) 



J^ < e 



(r a - 2 u 2 + r 2n - 2+a ul)dx 



':rt 
-M 



rlVl 

+C e (l + II {r n - l u) x \\l°°) / [r 2n - 2+a u 2 x + r a ~ 2 u 2 + r a {p - p c 
x Jo 



dx. (3.66) 



From (|3.12p - (|3.13p . using the same argument in the proof of (|3.44p and the Cauchy-Schwarz 
inequality, we get 



J2 + J3 





rM 



(2 C1 + c 2 )p 1+e (r n - l u t ) x {r n -' +a u t ) x - 2c x (n - l)p ti (r 



n— 1+a 



„ra-2+a„,2 



t Jx 



dx 



+ / {(2d + 02) (1 + e)p 9+2 [(r n - 1 u) x ] 2 - (n - l)(2d + c 2 )p 1+e ( 
Jo 

_ 7/0 7+i (r n-i n)x _ 2ci(n _ l)0 p 6+i( r n-l u)x ^ _ 2cx {n - 1)/^} 

dx 



1+e <r n ~ 2 u 2 ) 



P 



rM 

+2d(n-l) J |(n- l)r n ~ 2+a up 6 



r / x 



u t 
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2 

r J x \r 



< -(C 22 - e) / (r 2n - 2+a u^ + r a - 2 u 2 t )dx 
Jo 

r-M 

+C t (l + \\{r n - l u) x \\ 2 L¥ ) / [r 2n - 2 + a ul + r a ~ 2 u 2 + r a (p- PoD ) 2 ]dx, (3.67) 

•/ 



pM r-M 

J a < e 



r-IVl r-IVl 

/ r a ~ 2 u 2 t dx + C e ((l-n)Gxr- n u + d r Afu + d t Af) 2 r 2+a dx 
Jo Jo 

r-M pM 

/ r a - 2 u 2 dxds + C e (/| + / r 2+t Vdx) (3.68) 
io io 



and 



/•Af 

Js = - [(AP)'(r n - 1+a u t ) x + (n - ljA?^- 2 ^),] 

< e/ (r a - 2 u 2 + r 2n - 2+a u 2 xt )dx + C t (\AP\ 2 + \(APy\ 2 ). (3.69) 



Let e = ±C 22 , from (13T65J) - (13391) . we have 

-M /-M 



|jf r«u 2 t dx + ±fj Q (r 2 - 2 +«u 2 xt + r«- 2 u 2 t )dx 

r-M 

< C(l + {{(r^uUHoo) / [r 2n - 2+Q n 2 +r Q -V + r a (p- Poo ) 2 ] dx 
+C(/ 2 + |AP| 2 + KAP)'| 2 )+C 23ei ^ j( r- 2 u 2 dx. 
If C 23 ei < i from (I3T35|) and (13361) . we can obtain 

J r-M (~i r-M 

|jf r^dx + ^f^ (r 2 - 2 +«u 2 xt + r«- 2 u 2 t )dx 



r-lVl 

< C(l+[|(r n ~ 1 u)*ll£<») / [r^-^ + r^-V+r^p-^) 2 ]^ 
1 io 

+C(/ 2 2 + |AP| 2 + |(AP)'| 2 ) (3.70) 



and 



/ (r a u 2 )(x,t)dx + / / (r 2n - 2+Q u 2 4 + r Q - 2 u 2 ) (x, s)dxds 
io io io 

< Ce 2 (l + ||(r*- 1 u) a! ||£ s ). (3.71) 
From the equation (11.201) 9. we have 

(2 C1 + c 2 )r n ~ l d x {p l + e d x (r n - l u)) = u t + + 2c* (n - l)r n - 2 u(A + /, 

and using the estimates (pTTT]) . (^T2]) - (^T3| ). (|335l - (|337j ) and (f37?Tj) . conclude that 

M 

r 2n - 2+a [d x ( P 1+e d x (r n - l u))] 2 dx < Ce 2 (l + Mr^uUlg), 
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and 

cM 

\d x {p l+e d x {r n - l u))\dx < Ce (l + W^uU^), (3.72) 

for all t G [0,T]. From (]3TT2]) , (pT54"|) and (f3772]) . using Sobolev's embedding Theorem ^ 
we can obtain 

Wdx^-^hs < C 24 e (l + Kt^uUl-). (3.73) 

If C 24 e < 5, from (13341) . (|37TL]) and (pT73]l . we can get (pr62>(pO)4j) immediately. □ 
Now, we can let 

ei = (Ci6 + C 20 )e . (3.74) 

If (1 + | + g + g| + 4C 8 + 2Cf 4 + 4Ci 7 + 2C 23 )ei + 2C 24 e < 1, using the results in Lemmas 
I3.3M3.81 we finish the proof of the Claim 1. Thus, we can let eo be a positive constant satisfying 

(1 + - + ^ + ^ + AC 8 + 2Cf A + 4C 17 + 2C 23 )(C 16 + C 20 )e + 2C 24 e = 1. (3.75) 
Using the classical continuity method, we can obtain the following lemma. 

Lemma 3.9. Under the assumptions in Theorem the solution (p,u) satisfies the estimates 

Remark 3.1. In this section, we just give some ideas. In deed, we use the continuity method to 
obtain the bound of ^r^j, i = 0, . . . , N + 1, in Section 2J The basic theory of ordinary differential 
equations guarantees /?;(£), Uj(t), r;(t) G C([0, T*)), i = 0, . . . , N. Since rj(t) > h, j = 0, . . . , N, 
we have G C([0, T*)), z = 0, . . . , iV. Thus, we can use the continuity method. 
From Lemma ESI we can obtain the following lemma easily. 

Lemma 3.10. Under the assumptions in Theorem \l.ll if eo is small enough, we have 

\\p(;h) - p(;t 2 )\\ L 2 <C\t 1 -t 2 \, 

IK,ti) -u{-,t 2 )\\ L 2 < c\h-t 2 \, 

||r(.,ti)-r(.^ 2 )IUoo < C\h-t 2 \, 
\\d x {r n ~\)(-M)-d x (r n - 1 u)(-M)\\^<C\t l -t 2 \^, 



„n— 2„,\ ( „n—l\ 



c )(;t)\\ L n-^<C, 



for all ti,t 2 ,t > 0. 



4 Difference scheme and approximate solutions. 

In this section, applying a discrete difference scheme as in [TJ, we construct approximate solutions 
to the initial boundary value problem (|1.20|) - (|1.23|) . 

For any given positive integer N, let h = i be an increment in x and xj = jh for j G 
{0, . . . , N}. For each integer N, we construct the following time-dependent functions: 

(^(t),^),^)), j = 0,...,N, 

that form a discrete approximation to (p, u, r)(xj,t) for j = 0, . . . , N. 
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First, pi(t), Uj(t) and rj+i(i), i = 0, . . . , N, j = 1, . . . , N, are determined by the following 
system of 3iV + 2 differential equations: 

j fPi = -p^rr^il (4-1) 

± Uj = r^Saj - 2(n - l)r?" V(^-i) - fj, (4.2) 

4^+1 = Uj+i, (4.3) 

with the boundary conditions: 

uo(t) = 0, r$(t) = h, (4.4) 

Pn ~ Pn(Xn + 2p N )5(r n N 1 u N ) + 2(n - 1)^-/^ = ft, (4.5) 

?"jV+l 



and initial data 



(Pi,«i)(0)= f \ T Pv{y)dy,\ F u {y)dy ) ,j = 1,...,N, (4.6) 



h.l(j-i)h h J{j-i)h 



Po(0) = pi(0), «o(0) = 0, rff(0) = fc, (4.7) 
*-i , 

r-(0) = / l + n^— - i = l,...,iV + l, (4.8) 
«=o ^ ' 

and ii7v + i(0) satisfies 

ft(0) - Pn(0)(Xn(0) + 2 m (0))5(r"~ 1 (0)n w (0)) + 2(n - l) UN+l( ;% N (0) = ft(0), (4.9) 

rAT + i(0) 

where <5 is the operator defined by Swj = {wj + % — wj)/h, and 

aj{t) = /9j_i(Aj_i + 2 / u j _i)5(r™r 1 1 'u j _i) - Pj_i, 

A i = Kpj),H = P(pj), p j = P(Pj), 
f j (t) = f(jh,r j ,t). 

The boundary conditions (|4.4[) - (|4.5)) are consistent with the initial data. The condition (|4.5p 
determines UN+i(t). 

Let (poo.i^So^) = (Poo(ih), h + r^(ih)), i = 0, . . . , iV, we have 

^(^2oj-i) = -^ + Oy» 

i-i ^ 

fc=o Poo ' fc 

and 

1— n 2 

iQiil < C(jh)— h, \Q 2j \<C(jh)^h, j = l,...,N. 
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Then, for any small h, the initial data (j>i, Ui, r.j)(0) and the external force /j, i = 0, . . . , N, 
satisfies 



N-i 

^max n \(h(P) - PooM 2 + E [ r T ~ 2+a ( S Pi ~ S Poc,j) 2 } (0)h < Ce 2 , 



ie{o,...,JV} 



(4.10) 



3=0 



N 

C-\i + l)h < rf (0) < C(t + E [r7\| + rf - 2 (^) 2 ] (0)/i < Ce 2 ,, (4.11) 



i=o 



N 

E 



,2n-2a 



M^(^-iS-i)) )(o)^<c4 



(4.12) 



where C > are independent of /i. 

The basic theory of differential equations guarantees the local existence of smooth solutions 
(pi, Ui, Ti) (i = 0, . . . , N) to the Cauchy problem (|4.ip - (|4.9p on an interval [0, T h ), such that 

< pi(t) < do, |n(i)| < oo,i = 0,... ,JV, 

with the aid of (l4TTU|) - (|4TT2l . 

For any fixed T > 0, by virtue of Lemma 13. 1113.101 and using similar arguments as in [HE], 
we can obtain the following lemma and prove that the Cauchy problem (|4.ip - (|4.9p has a unique 
solution for t G [0, T] when h < hTe i where fyre > is a constant dependent on T and eo- 

Lemma 4.1. For any h G (0, fyr,e L there exist a positive constant C independent of h such that 



Pi(t) G 



1 3_ 
2* 2 P 



\pi{t) — Poo,i\\l°° < Ceo, 



r?(i) G [C-\l + l)h,C{l + l)hl 



N 



Y,(u 2 3(t) + \Plit)-Pooj\ 2 )h<Cel 

3=0 

Ul{t) 



n{t) 



<Ce , 



N 



E K + »f ~ 2 (H) 2 + 4 < Cel 

j=0 \ T 3 ) 

e (I + -r 2 (^) 2 ) + jf e (^) 2 (s)Ms ^ Ce 



/ E i r j(P3 - Poo,jf + (rj ~ roojf] hds < Ce 
Jo J=0 



v 



3=0 



N 



^(rfv$)(t)dx + I ( rf-\8 Uj f + -L}hds< Ce^ 



3=0 



3 



N-1 ft N-1 

E (»f"* fa (*Pi - + / E (*f" a+a (*Pi - W)WW» < Ce o 

i=o - 70 j=0 
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N-l 

\Spj - 5pooj\(t)h < Ce , 

3=0 

N r / a \ 2 ~ 

3=0 L V 7 . 



E i\pj(h) ~ Pj{t2)\ 2 + \uj(h) - u 3 {t 2 )\ 2 ) h < C\h - t 2 \ 2 , 

3=0 

\n(h) - n(t 2 )\ < c\h-t 2 \, 



N 



\5{r]- l Ui){t x ) - S^uAfotfh < C\h - t 2 \, 



3=0 

N 

E (i^rsor 1 + i^rr 1 )!""") h < c > 

3=0 

for alltx,t 2 ,te [0,T], ie{0,...,N} and I G {1,...,JV + 1}. 

Now, we can define our approximate solutions (p N ,u N ,r N )(x,t) for the Cauchy problem 
(jl.20p - ([1.23p . For each fixed ./V and t G [0, T], we define piecewise linear continuous functions 
(p N ,u N ,r N )(x,t) with respect to x as follows: when x G [[#-ZV], [xN] + 1] 

p N (x,t) = p [xN] (t) + (xN - [xN])(p [xN]+1 (t) - p [xN] (t)), 
u N {x,t) = u [xN] (t) + (xN - [xN])(u [xN]+1 (t) -u [xN] (t)), 

r N (*,t) = (rf xN] (t) + (xN - [xN])(r? xN]+1 (t) - if aJV] (0)) V " . 

From Lemma 14. 1\ using similar arguments as in [BE], we can obtain the compactness of 
approximate solutions (p N ,u N ,r N ) and prove the existence part of Theorem 11.11 Since the 
constant C in Lemma 14.11 is independent of T, we can obtain the regularity estimates (|1.29j) - 
(fOTp easily. 



5 Uniqueness. 

In this section, applying energy method, we will prove the uniqueness of the solution in Theorem 
11.11 Let (pi,ui,ri)(x,t) and (p 2 ,u 2 ,r 2 )(x,t) be two solutions in Theorem 11.11 Then, we have, 
i = 1,2, (x,t) G [0,M] x [0,2*1 

, C~ x x^ < n(x,t) < Cxn, (5.1) 

\x~"Ui(x,t)\ + \x~~n~ d x Ui(x,t)\ < C. (5-2) 

For simplicity, we may assume that (pi,ui,ri)(x,t) and (p 2 ,u 2 ,r 2 )(x,t) are suitably smooth 
since the following estimates are valid for the solutions with the regularity indicated in Theorem 
11.11 by using the Friedrichs mollifier. 



Pi(x,t) G 



2* 2 P 
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Let 

Prom (|3.ip . we have 



q = p\ — p 2 , w = u\ — U2, R = r± — r 2 . 



d 
dt 



f M -1 2, \ f M -1 

/ x ™ R (x,t)dt = 2 / x n Rwdx 
Jo Jo 



< e 



M 



X n W d,X + C € 



M 



— — 1 

x nR dx. 



From (fL20j) and (|ET ]) - ([Q]) . we have 



dt 
2 



M 



M 



q (x,t)dt = 2 gdt(pi — P2)dx 



M 



o | -pfr^d^ + p\rr x d x u 2 - (n - 1)^ + (n - 1) 

T\ r 2 



P2U2 



f 2n~2 2 f 2 

< e / (x~ ~ w 2 + x~~w 2 )dx + C 6 / (g 2 + x~™ R 2 )dx. 
Jo Jo 

From the equation ()1.2Qj) ? and boundary conditions (|1.22[) - (|1.23j) . we get 
d /* M 1 2 ( ^ J 

dtj 2 W{x > t)dx 



At 



+ {{-c l + C2)p\ +6 [(r n 1 - 1 w) x ] 2 + 



+6\un-l„„\ 12 , 2(n — 1) 1+0 



n 



■n 



a^r™- 1 ^) (2d + c 2 )(p\+ e - pl+^rr 1 ^) 



+(2 C1 + c 2 )^ +e ^((rr 1 - rr 1 )^) - (PT - pi) 

r M r 111 

+ / 2d(n-l)<9 :E r^uW ) p\dx 

Jo L r i r 2 J 

r M r 1 

+ j( 2d(n - 1)^ rT x w{f) (p? - ^)dx 



r/.r 



/■M 

+ / 2 Cl (n-l)p^ 

JO 



n-1 ,„n-:h„. n2 



''2 



r )w- 



dx 



M 



ASr((rJ 



n— 1 _ra— 2\ 



(2c 1 +c 2 )pi + ^(rr 1 n 2 )-pl 



c/.r 



+ 



wGx(r 



l—n „1— n 



1 1 - n )dx+ / W (A/(x,r 2 ,t)-A/( 2 ;,r 1 ,t))dx. 
JO 



From (EB-ES) and (1531) . we have 



ri 



fit 

< C 



M 



1 



^ / ^w 2 (x,t)cix + C 22 



2n^2 2 _2 2 -l 

x » w^ + i n w > ax 



M 







x ~R 2 + g 2 + w 2 ) dx. 



From and (j5U|) . letting e = ±C 22 , we obtain 

dt 



r M 2 r } 

\ [w 2 + g 2 + X -nR 2 ]dx < C 
Jo Jo 



x ™R 2 + q 2 + w 2 ) dx. 
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Using Gronwall's inequality, we have for any t G [0, T], 



f M 2 
/ [w 2 + g 2 + x~nR 2 ]dx = 0. 
Jo 



This prove the uniqueness of solution in Theorem 11,1 



6 Asymptotic behavior 

In this section, we consider the asymptotic behavior of the solution to the free boundary problem 
(ll,20p - (jl.23p . We will show that the solution to the free boundary problem tends to the stationary 
solution as t — » +oo. 

The following lemma is proved in |15j . 

Lemma 6.1. Suppose that y G W^(M. + ) satisfies 



and 



V = 2/1 + 2/2, 



1 2/2 1 < J2 ai > - 



on 



i=l 



i=l 



where y\ G PF, ' _(R + ), and lim 2/1 (s) = and aj,/3j G L Pi 

1, ...,m. Then lim y(s) = 0. 

s— >+oc 



for some pi G [1, oo), i 



Proposition 6.1. Under the assumptions of Theorem \l.l\ the total kinetic energy 



f M 1 

E(t) := / -u 2 (x,t)dx — > as t — >■ +oo. 
io 2 



Proof. From f)3. 15j) and Lemma [3.9l we have E(t) G L 1 (R + ). Using the Cauchy-Schwarz inequal- 
ity, we obtain 



\E'(t)\ < E(t) + / u 2 t dx. 
Jo 



Taking into account the estimate (|3.64|) and Lemma 13.91 applying Lemma 16. 1\ we finish the 
proof. □ 

Proposition 6.2. Under the assumptions of Theorem \l.l[ we have 

i-M 



I-Nl 

I {t — r OD ) 2 (x,t)dx — > as t — > +oo. 
Jo 



Proof. From (|3.20p and Lemma 13.9^ we have J (r — r 00 ) 2 (x,t)dx G L 1 (M + ). Using the Cauchy- 
Schwarz inequality, we obtain 



d f M 

37 / (r - rcofdx 
dt Jo 



< 2E(t) + / (r - roD ) 2 dx. 
Jo 

Taking into account the estimate E(t) G L 1 (M + ), applying Lemma 16. 11 we finish the proof. □ 



2 / (r — r^udx 
o 
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Proposition 6.3. Under the assumptions of Theorem \l.l[ we have 

-M 



and 



as t —* +00. 



(p ~ Poo) 2 (x,t)dx -> 0, 
\\(p- Poo)(-,t)\\L« ->0, «G [l,oo), 



(6.1) 
(6.2) 



Proof. From (|3.2U|) and Lemma 13.91 we have J (p — p OQ ) 2 (x,t)dx G L X (]R + ). From (jl.29|) . using 
the Cauchy-Schwarz inequality, we obtain 



d_ 
dt 



M 



(p - Poofdx 



M 



< (P~ PooYdx + C / (r n - A u)^da:. 



Taking into account the estimate (|3.15p and Lemma 13.91 applying Lemma l6.lt we obtain (|6.1|) . 
From (fL29l) . (1331) and (|53]> . we can obtain flO]) easily. □ 

Proposition 6.4. Under the assumptions of Theorem \l.l\ we have 



M 



x n {(p )x-(Poc)x) (x,t)dx^0, as t -> +00. 







Proof. From (fL29l . (f335D . (|373T1) and Lemma we have 



x =! ^(( P ) x - (p^fix^dx G L 1 ^)- 







From (|1.20p . (|1.29p - (|1.30|) and f|3.3[) . using the Cauchy-Schwarz inequality, we have 



d 
dt 

20 



x -ir-((p») x -( P l) x fdx 



M 



2n-2 + a 



20 



x^^((p% - (pl) x )(p» +1 d x (r n -t U )) x dx 

M 

x^({p\-(pl) x ) (^ + A(p% 







2ci + c 2 

,u(p) x f(x,r,t) , , 
+2ci(n - 1) VP 1 + Jy ] \ ' ) dx 



r n—l 



< c 



M r 



«P°)» ~ (PL),) 2 + r a ui + r a (r - r^) 



+r a (p - Poo y + r Q " VJ dx + ft 

Taking into account the estimates (jl.29p . (|3.35p - (|3.37p . (|3.62[) and Lemma [3791 applying Lemma 
6-H we end the proof. □ 



From Proposition I6.3lf6.4l using Sobolev's embedding Theorem, we can obtain the following 
corollary immediately. 

Corollary 6.1. Under the assumptions of Theorem \l.l\ we have 

\\p(-,t) - P oo(-)\\l°° + \\r(-,t) - roo(-)||i/x> 0, ast^ +00. 
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Proposition 6.5. Under the assumptions of Theorem \l.l[ we have 

2ra — 2+g q / 

x n u x (x, t)dx — ► 0, (IS t — ► +00. 



Proof. From the estimates (jl.29j) . ()3.36p and Lemma 13, 9\ we have 



M 



x 2n n +a u 2 x {x,t)dx G L 1 (R + ) 



Using the Cauchy-Schwarz inequality, we have 



dt 



M 



2n — 2 + q o 
X n U„.dx 



M 



2n-2 + q 

X « U x U xt dx 







< 



A7 



Ai 



. 2n-2+q 3 

u^cta + / x n u x tdx. 
'0 



Taking into account the estimates (|1.29p , (|3.36p , (|3.62p and Lemma 13.91 applying Lemma 16.11 
we end the proof. □ 

From Proposition 16.11 and 16.51 using Sobolev's embedding Theorem, we can obtain the fol- 
lowing corollary immediately. 

Corollary 6.2. Under the assumptions of Theorem \l.l\ we have 

||tt(-,t)||ioo — ► 0, as t — > +oo. 
Thus, we finish the proof of Theorem 11.11 



7 Stabilization rate estimates 

Now we are in position to estimate the stabilization rate. We first state the following proposition 
which gives the stabilization rate estimates in L 2 ([0, M])-norm of the solution. 



pM 

/ (n 2 + {p- Poo f + x~\r n - r£) 2 ) dx < Ce 
Jo 



ait 



Proposition 7.1. Under the assumptions of Theorem \1.3l we have 

10 

and 

\p(M,t)- Poo (M)\+n r' n --\p- Poo )ldx\ + \\r(.,t)-r OQ (x)\\ L2 <Ce-^ 

for all t > 0, where a\ is a positive constant. 
Proof. Let 

r M i 

V x = / -u 2 dx + S[V}- SlVov], 



ait 



W 1 = J^r^u)l + r^ul + ^dx. 



(7.1) 



(7.2) 
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From (Q2l) . ([3TT8>([3TT8l) . we have 



V{ + 2C 31 iyi < CfiV{ 2 + C7|AP| 2 < Ce~ aot V~ 2 + Ce 



2aot 



M 



Cm / (« 2 + (P " Poo) 2 + *~V - r^f) dx 



and 

From (pT28j) . we have 



< Vi < C 32 / (n 2 + (p - Poo ) 2 + x- V " ^) 2 ) ^, 

JO 



- 2 (r"-r^) 2 l dx 



/ [(p-poo) 2 +x- 2 (r 
JO 

_ ( !^ _ 1°° ) dx + CssWr + Ce~ 2aot . 



~ - C38 Jt 



d f M u I ' r n r" 



n n 



From (|1.29p . we obtain 

C38 

Let 



' ' 00 1 dx 



r n 1 V " " 



^ dx < C39 J (u 2 + \p - poo\ 2 ) dx. 



f M v fr n r n 

V 2 = V 1 + eC 38 20)dx, 

Jo r \ n n 



W 2 = C 31 Wt + e / [(pT - p2o) 2 + x" V - C) 2 ] dx, 
Jo 



where e = mmj&J-, 2C 1 C }. From (|T.3j) and (|7.4p - (|7.7p . we have 



V^ + W 2 < Ce~ 2aot , 



M 



C39 1 / («' + (P " Poo) 2 + *~ V - r^f) dx 



< 



M 



V 2 < C39 / (n 2 + (p - poo) 2 + x~V " O 2 ) dx, 
Jo 



and 



C 40 / (« 2 + (P - Poo) 2 + x~ 2 {r n - r^) 2 ) dx < W 2 . 
Jo 



(7.3) 



(7.4) 
(7.5) 



(7.6) 



(7.7) 



(7.8) 



(7.9) 



(7.10) 



Thus V 2 is a Lyapunov functional. From (jl.32p . we obtain the estimate (|7. 1|) . From (jl.29p . 
(j535D , ([530D and dH}, we can get ([721) easil y- □ 



Proposition 7.2. Under the assumptions of Theorem \1.3l we obtain 



M /y 2 



2n—2 2 
_|_ r *n ^ 



10 \r* 

for all t > 0, where 03 is a positive constant. 



2 }j (x,t)dx < Ce~ a3t , 



(7.11) 
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Proof. Let 



2 



J [2 n n rp 

MApl, - Ap^ + &P)(r n - l u) x + u (g-^ - G^^j \ dx. 

From (fPOj) and (laMI) . we have 

r M 1 

Vl + J -u 2 (x,s)dx<C il {f 2 + \(APy\ 2 + W 2 ). 

From (|3TL3l) - ([3TT4"1) . we have 

rM 2 
V 3 > / {C, 2 {(r n - l uf x + ^ + r 2n ~ 2 u 2 x ) - C 43 ((p - Poo ) 2 + |AP| 2 )}^, 
Jo r 

r M ii 2 

V 3 < / {C^({r n ~ l u)l + - + r 2n - 2 u 2 x ) + C 43 ((p - Poo) 2 + |AP| 2 )}dx. 
Letting V4 = V 2 + t]V 3 + r)C 43 \AP\ 2 , where r, = min{±, 4C3 * C43 , Fr om iTT^- imPl) . we have 

CW 2 >V 4 > C~ l W x , 

and 

y 4 ' + C- X W 2 < C (f 2 + | AP| 2 + |(AP)'| 2 ) . 
Thus V4 is a Lyapunov functional. From (jl.32p . we can obtain the estimate (|7.1ip , □ 

Proposition 7.3. Under the assumptions of Theorem \1.3l we obtain 

rM , 

Jo 



r2~ m (p - Poo ydxds < Ce~ at , (7.12) 



M 

rh- m ( r - r^fdxds < Ce~ at , (7.13) 







AI 



1 



(r2- m u 2 )(x,t)dx < Ce~ at , (7.14) 







M 

(r 2n - 2+ 2- m (p - Poo ) 2 x ){x, t)dx < Ce~ at , (7.15) 

for all t > and m = 0, 1, . . . , n — 1, where a is a positive constant. 

Proof. From (|7.1[) - (|7.2|) . we know that the estimates (|7.12|) - (|7.15|) hold with m = 0. 

Claim 3: If that (f77T2]) - (f77T5|) hold with m<k,k<E [0, n-2], then the estimates ([77T2]) - (f77T5]) 
hold with m = k + 1. 

We could prove Claim 3 as follows. Let = ^ — k — 1. From (|3.4ip . (|7.ip and (|7.15p (m = A;), 
we have 

rM 

/ r^(/>-Poo) 2 ^<Ce- at , 
Jo 

and (f7TT2|) (m = k + 1) holds. From (fl~29|) and ([7TT2]) . we can obtain ([7TT3]) (m = fc + 1) easily. 
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From (|3.49p . we obtain 



d_ 

di 



,-M r M / 2\ 

J r° k u 2 dx + C 42 J r a * \ r 2n - 2 ul + ^jdx<C (ff + \AP\ 2 ) + Ce""'. (7.16) 



Thus /q (r ak u 2 )(x,t)dx is a Lyapunov functional, and we obtain (|7.14p (m = fe+1) immediately. 
From (13321) and (ffl^ - fflil) . we have 



d_ 

dt 



/ H 2 dx + ^ / H 2 cte 
/•M 

< C / (r a «u 2 + r a <=(p - p^) 2 + r afc (r - ^co) 2 ) + Cff < CeT 01 
Jo 



Thus Jq M H 2 (x,t)dx is a Lyapunov functional. Using the estimates (jl.29p and (|7.12p - (|7.14p . we 
obtain (|7.15p (m = k + 1), finish the proof of Claim 3 and Proposition 17.31 immediately. □ 

From (jl.29p . (|7.12p and (|7.15p . using Holder's inequality and Sobolev's embedding Theorem, 
we could obtain the following proposition. 

Proposition 7.4. Under the assumptions of Theorem \1.3l we obtain 

\\p(;t) - PooQIU- + \H;t) - r^-)]^ < Ce~ at , 
for all t >0, where a is a positive constant. 

Proposition 7.5. Under the assumptions of Theorem \1.S\ we obtain 

M 

(r a u 2 )(x,t)dx <Ce~ a \ (7.17) 



o 

for all t > 0, where a = | — n and a is a positive constant. 
Proof. From (jl.30p and (|3.70p . we have 



| J r a u 2 dx + C i5 J o (r 2n - 2+a u 2 xt + r^ 2 u 2 ) dx 
pM 

< C m [r 2n - 2+a u 2 x + r a - 2 u 2 +r a (p- Poo ) 2 }dx 
Jo 



+C (fl + \AP\ 2 + \{AP)'\ 2 ) . (7.18) 

Let V A = J^(r a u 2 )(x,t)dx + J Q M (r a u 2 )(x,t)dx. From (fL32j) . (ITTZI) . (|71Sj> (ifc = n - 2) and 
dSSD , we have 

F 4 ' + C^Va < Ce~ aot + Ce~ at . 
Thus ^ is a Lyapunov functional, and we obtain (17.171) immediately. □ 
Proposition 7.6. Under the assumptions of Theorem 1 1 . 31 w;e obtain 



J r 2 "- 2+Q [9 :c (p 1+9 ^(r"- 1 U ))] 2 dx+ K^K" 1 ^) (■,*) 
/or a// i > ; where a = \ — n and a is a positive constant. 



< Ce~ at , 

L°° 
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Proof. From the equation (jl.20p ? . we have 

(2 C1 + c 2 )r n - 1 d x (p l+e d x (r n - l u)) = u t + Ar"- V)* + 2c i( n ~ l)r n ~ 2 u(P% + f, 
and using the estimates (O9l) - (fL30j ). (Q2"1) . (l7TT2ll - (l7TT5D and (I7TT71) . conclude that 

r 2n - 2+a [d x (p 1+e d x (r n - 1 u))] 2 dx < Ce~ at , 



o 



and 

pM 

/ |^(p 1+ ^(r n ~M)l^ < Ce- at . (7.19) 
J o 

From (jl.29p . (|7.1ip and (|7.19p . using Sobolev's embedding Theorem W 1 ' 1 L°°, we can obtain 

\\d x (r»- x u)(;t)\\L~> <Ce~ at , 

and 



U / X 

- •,*) 
r 



L°° 



□ 
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